Third Newton Law in the Mean and Perturbations of Noise 
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We analyze an isothermal Brownian motion of particle 
ensembles in the Smoluchowski approximation. Contrary to 
standard procedures, the environmental recoil effects associ- 
ated with locally induced heat flows are not completely disre- 
garded. The main technical tool in the study of such weakly- 
out-of equilibrium systems is a consequent exploitation of the 
Hamilton- Jacobi equation. The third Newton law in the mean 
is utilised to generate diffusion-type processes which are either 
anomalous (enhanced), or generically non-dispersive. 



If we consider a fluid in thermal equilibrium as the 
noise carrier, a kinetic theory viewpoint amounts to vi- 
sualizing the constituent molecules that collide not only 
with each other but also with the tagged (colloidal) par- 
ticle, so enforcing and maintaining its observed erratic 
motion. The Smoluchowski approximation takes us away 
from those kinetic theory intuitions by projecting the 
phase-space theory of random motions into its configura- 
tion space image which is a spatial Markovian diffusion 
process: 
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In the above m stands for the mass of a diffusing par- 
ticle, (3 is a friction parameter, D is a diffusion constant 
and W{t) is a normalised Wiener process. The Smolu- 
chowski forward drift can be traced back to a presumed 
selective action of the external force F — — VU on the 
Brownian particle that has a negligible effect on the ther- 
mal bath but in view of frictional resistance imparts to a 
particle the mean velocity F/m/3 on the /3 _1 time scale, 
@§ 

Smoluchowski diffusions are conventionally regarded as 
isothermal processes. If we however admit that tiny heat 
flows should accompany the particle transport due to the 
Brownian motion, a suitable description of thermal inho- 
mogeneities and their effects on the dispersion of Brow- 
nian particles must be invented. 
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To exemplify this point, let us consider the phenomeno- 
logical assumption of Ref. || : "a gas of Brownian parti- 
cles falling in gravity should leave a trail of warm fluid 
in its wake, since its potential energy is being converted 
into heat" . Obviously, if the particles would move against 
gravitational force, then the temperature of the medium 
should locally drop down. Those features, if we are to 
keep track of the local heating and cooling (as opposed 
to the isothermal Einstein or Smoluchowski diffusive dy- 
namics) were interpreted in Ref. || as a source of the 
space-time dependence of temperature. In that case the 
Fokker-Planck equation must be supplemented by an evo- 
lution equation for the temperature field. 

The phase-space scenario usually refers to minute ac- 
celeration/deceleration events which modify (say, at a 
rate of 10 21 per second) velocities of realistic particles. 
Clearly, the microscopic energy-momentum conservation 
laws need to be respected in each separate collision event. 
In contrast to derivations based on the Boltzmann col- 
llision scenario, this feature is completely alien to the 
Brownian motion theory, cf. [™q|. Therefore, there seems 
tempting to require that each minute acceleration of a 
Brownian particle is accompanied by a minute cooling 
of the medium in its immediate neighbourhood, while 
any deceleration event should induce a local heating phe- 
nomenon, see e.g. Q]. Since (cf. Eq. (1)), we always 
disregard the fine details of about 10 13 collision impacts 
on the Brownian particle on a typical relaxation time 
scale of 10" 8 s, there is definitely enough room to allow 
for local statistical measures of heating and cooling. 

In a weakly out-of-equilibrium system with small heat 
flows, we may expect that the standard equilibrium tem- 
perature notion needs a replacement by an effective tem- 
perature notion (and an effective thermal equilibrium), 
which depends on the chosen fast-versus-slow-process 
time scales and the ensemble averaging. In the above 
sense only an effective isothermal regime may be main- 
tained. 

It is well known that a spatial diffusion (Smoluchowski) 
approximation of the phase-space process, allows to re- 
duce the number of independent local conservation laws 
(cf . |^-^| ) to two only. Therefore the Fokker-Planck equa- 
tion can always be supplemented by another (indepen- 
dent) partial differential equation to form a closed sys- 
tem. Non-isothermal flow description needs to accomo- 
date variations of temperature of the bath, (cf. ||), while 
we investigate the limits of validity of the isothermal sce- 
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nario. That amounts to inequivalent choices of the sup- 
plementary equation. 

If we assign a probability density po(x) with which 
the initial data xq = X(0) for Eq. (1) are distributed, 
then the emergent Fick law would reveal a statistical ten- 
dency of particles to flow away from higher probability 
residence areas. This feature is encoded in the corre- 
sponding Fokker-Planck equation: 

d tP = -V • (vp) = -V • [(-^ - D^)p] (2) 

mp p 

where a diffusion current velocity is v(x, t) = b(x, t) — 
D Vp S-?y while the forward drift reads b(x,t) = 
Clearly, the local diffusion current (a local flow that 
might be experimentally observed for a cloud of sus- 
pended particles in a liquid) j — vp is nonzero in 
the nonequilibrium situation and a non-negligible mat- 
ter transport occurs as a consequence of the Brownian 
motion, on the ensemble average. 

It is interesting to notice that the local velocity field 
v(x, t) obeys the natural (local) momentum conservation 
law which directly originates from the rules of the Ito 
calculus for Markovian diffusion processes, [pj, and from 
the first moment equation in the diffusion approximation 
(!) of the Kramers theory, ||: 

d t v + (v- V)i? = V(0 - Q) . (3) 

An effective potential function f2(x) can be expressed 
in terms of the forward drift b(x) = — ^ as follows: £1 = 

F 2 + _2_V • F 

2m 2 /3 2 ^ m/3 V r • 

Let us emphasize that it is the diffusion (Smolu- 
chowski) approximation which makes the right-hand-side 
of Eq. (3) substantially different from the usual moment 
equations appropriate for the Brownian motion. In par- 
ticular, the force F presumed to act upon an individual 
particle, does not give rise in Eq. (3) to the expression 
— — Vy which might be expected on the basis of kinetic 
theory intuitions and moment identities directly deriv- 
able from the Karmers equation, but to the term 

Moreover, instead of the standard pressure term, there 
appears a contribution from a probability density p- 
dependent potential Q(x,t). It is given in terms of the 
so-called osmotic velocity field u(x,t) = DV lnp(x,t), 
(cf. [H): Q(x, t) = ^u 2 + DV ■ u and is generic to a lo- 
cal momentum conservation law respected by isothermal 
Markovian diffusion processes, cf. |^,|| 

To analyze perturbations (and flows) of the noise car- 
rier, a more general function b(X(t),t), must replace the 
Smoluchowski drift of Eqs. (1), (2). The forward drifts 
modify additively the pure noise (Wiener process entry) 
term in the Ito equations. Under suitable restrictions, we 
can relate probability measures corresponding to differ- 
ent (in terms of forward drifts !) Fokker-Planck equations 
and processes by means of the Cameron-Martin-Girsanov 



theory of measure transformations. The Radon-Nikodym 
derivative of measures is here involved and for suitable 
forward drifts which are are gradient fields that yields, [||, 
the most general form of an auxiliary potential f2(x, t) in 
Eq. (3): 

n(a?,t) = 22W+i(^ + V-?)] . (4) 

We denote b(x, t) = 2DV(j)(x, t). 

Eq. (4) is a trivial identity, if we take for granted that 
all drifts are known from the beginning, like in case of 
typical Smoluchowski diffusions where the external force 
F is a priori postulated. We can proceed otherwise and, 
on the contrary, one can depart from a suitably chosen 
space-time dependent function fl(x, t). From this point of 
view, while developing the formalism, one should decide 
what is a quantity of a primary physical interest: the 
field of drifts b(x,t) or the potential f2(x, t). 

Mathematical features of the formalism appear to de- 
pend crucially on the properties (like continuity, local 
and global boundedness, Rellich class) of the potential 
f2, see e.g. ||[)|. Let us consider a bounded from below 
(local boundedness from above is useful as well), contin- 
uous function Q(x,t), cf. Q). Then, by means of the 
gradient field ansatz for the diffusion current velocity 
(v = V5 — > dtp = —V ■ [(VS')p]) we can transform the 
momentum conservation law (3) of a Markovian diffusion 
process to the universal Hamilton-Jacobi form: 

il = d t S+^\VS\ 2 + Q (5) 

where Q(x,t) was defined before. By applying the gra- 
dient operation to Eq. (5) we recover (3). In the above, 
the contribution due to Q is a direct consequence of an 
initial probability measure choice for the diffusion pro- 
cess, while il via Eq. (4) does account for an appropriate 
forward drift of the process. 

Thus, in the context of Markovian diffusion processes, 
we can consider a closed system of partial differential 
equations which comprises the continuity equation dtp = 
— V(t/p) and the Hamilton-Jacobi equation (5), plus suit- 
able initial (and/or boundary) data. The underlying dif- 
fusion process is specified uniquely, cf. 

Since the pertinent nonlinearly coupled system equa- 
tions looks discouraging, it is useful to mention that a lin- 
earisation of this problem is provided by a time-adjoint 
pair of generalised diffusion equations in the framework of 
the so-called Schrodinger boundary data problem, HQ. 
The standard heat equation appears as a very special 
case in this formalism. 

The local conservation law (3) acquires a direct phys- 
ical meaning (the rate of change of momentum carried 
by a locally co- moving with the flow volume, |6|), only if 
averaged with respect to p(x, t) over a simply connected 
spatial area. If V stands for a volume enclosed by a 
two-dimensional outward oriented surface dV, we define 
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a co-moving volume on small time scales, by deforming 
the boundary surface in accordance with the local cur- 
rent velocity field values. Then , let us consider at time t 
the displacement of the boundary surface dV{t) defined 
as follows: x £ dV — ► x + v(x,t)At for all x G dV . Up 
to the first order in At this guarantees the conservation 
of mass (probability measure) contained in V at time t i. 

e - fv(t+At) Pi* > 1 + At ) d3x - Sv(t) t ) d3x ~ °- 

The corresponding (to the leading order in At) quan- 
titative momentum rate-of-change measure reads, cf. ||, 

J v pV{n - Q)d 3 x. 

For a particular case of the free Brownian expansion 
of an initially given p Q (x) = {7r J )a/2 exp(-^ I ), where 

a 2 = 4Z?i , we would have — J v p\/Qd 3 x = — J gv Pda, 
where Q(x,t) — gnj^p — 2(t+t ) ' w ^ ne the "osmotic 
pressure" contribution reads P(x, t) = — 9 ^ tg ^ p(x, t) for 
all x G R 3 and t > 0. 

The current velocity v(x,t) = VS(x,t) = 2 {t+t ) IS 
linked to the Hamilton-Jacobi equation dtS + ^\\7S\ 2 + 

-2 

Q = whose solution reads: S(x,t) — 4 ( t + to ) + 
^Dln[4irD(t + 1 )]. 

Let us observe that the initial data vq = —D\7lnpo 

> i for the current velocity field indicate that we have 
totally ignored a crucial preliminary stage of the dynam- 
ics on the time scale, when the Brownian expansion 
of an initially static ensemble has been ignited and so 
particles have been ultimately set in motion. 

Notice also that our "osmotic expansion pressure" 
P(x, t) is not positive definite, in contrast to the familiar 
kinetic theory (equation of state) expression for the pres- 
sure P(x) = apP(x), a > appropriate for gases. The 
admissibility of the negative sign of the " pressure" func- 
tion encodes the fact that the Brownian evolving concen- 
tration of particles generically decompresses (blows up), 
instead of being compressed by the surrounding medium. 

The loss (in view of the "osmotic" migration) of mo- 
mentum stored in a control volume at a given time, 
may be here interpreted in terms of an acceleration 
— f v pX?Qd 3 x induced by a fictituous "attractive force". 

By invoking an explicit Hamilton-Jacobi connection 
(5), we may attribute to a diffusing Brownian ensemble 
the mean kinetic energy per unit of mass J v p\v 2 d 3 x. In 
view of < x 2 >= 6D(t + to), we have also f R3 p\v 2 d 3 x = 
-j^H^y. Notice that the mean energy J v pi^v 2 + Q)d 3 x 
needs not to be positive. Indeed, this expression iden- 
tically vanishes after extending integrations from V to 
R 3 . On the other hand the kinetic contribution, initially 
equal J R3 \pv 2 d 3 x = 3D/a 2 and evidently coming from 
nowhere, continually diminishes and is bound to disap- 
pear in the asymptotic t — > oo limit, when Brownian 
particles become uniformly distributed in space. 

Normally, diffusion processes yielding a nontrivial mat- 
ter transport (diffusion currents) are observed for a non- 
uniform concentration of colloidal particles which are re- 



garded as independent (non- interacting) . We can how- 
ever devise a thought (numerical) experiment that gives 
rise to a corresponding transport in terms of an ensemble 
of sample (and thus independent) Brownian motion real- 
isations on a fixed finite time interval, instead of consid- 
ering a multitude of them (migrating swarm of Brownian 
particles) simultaneously. 

Let us assume that "an effort" (hence, an energy loss) 
of the random medium, on the scale, to produce 
a local Brownian diffusion current [pv)(x,to) out of the 
initially static ensemble and thus to decompress (lower 
the blow-up tendency) an initial non-uniform probability 
distribution, results in the effective osmotic reaction of 
the random medium. This is the Brownian recoil effect 
ofRef. @. 

In that case, the particle swarm propagation scenario 
becomes entirely different from the standard one . First 
of all, the nonvanishing forward drift b — u is generated 
as a dynamical (effective, statistical here !) response of 
the bath to the enforced by the bath particle transport 
with the local velocity v = —u. Second, we need to ac- 
count for a parellel inversion of the pressure effects (com- 
pression +VQ should replace the decompression — VQ) 
in the respective local momentum conservation law. 

Those features can be secured through an explicit real- 
ization of the action-reaction principle which we promote 
to the status of the third Newton law in the mean. 

On the level of Eq. (3), once averaged over a fi- 
nite volume, we interpret the momentum per unit of 
mass rate-of-change J v pV(Cl — Q)d 3 x which occurs ex- 
clusively due to the Brownian expansion, to generate a 
counterbalancing rate-of-change tendency in the random 
medium. To account for the emerging forward drift and 
an obvious modification of the subsequent dynamics of 
an ensemble of particles, we re-define Eq. (3) by set- 
ting — J v pV(£l — Q)d 3 x in its right-hand-side instead 

of + J v pV(£l — Q)d 3 x . That amounts to an instan- 
taneous implementation of the third Newton law in the 
mean (action-reaction principle) in Eq. (3). 

Hence, the momentum conservation law for the pro- 
cess with a recoil (where the reaction term replaces the 
decompressive "action" term) would read: 

d t v + (v- V)v = V(Q-fi) (6) 

so that 

d t s + ^\vs\ 2 -Q = -n (7) 

stands for the corresponding Hamilton-Jacobi equation, 
cf. [7|JfC[|, instead of Eq. (5). A suitable adjustment (re- 
setting) of the initial data is here necessary. 

In the coarse-grained picture of motion we shall deal 
with a sequence of repeatable scenarios realised on the 
Smoluchowski process time scale At: the Brownian 
swarm expansion build-up is accompanied by the parallel 
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counterflow build-up, which in turn modifies the subse- 
quent stage of the Brownian swarm migration (being in- 
terpreted to modify the forward drift of the process) and 
the corresponding built-up anew counterflow. 

The new closed system of partial differential equa- 
tions refers to Markovian diffusion-type processes again, 
■ The link is particularly obvious if we observe that 
the new Hamilton- Jacobi equation (7) can be formally 
rewritten in the previous form (5) by introducing: 

n r = d t S + ^\VS\ 2 + Q (8) 

where Q r = 2Q — and represents the previously 
defined potential function of any Smoluchowski (or more 
general) diffusion process. 

It is Q r which via Eq. (4) would determine forward 
drifts of the Markovian diffusion process with a recoil. 
They must obey the Cameron-Martin-Girsanov identity 

fl r = 2Q - Q = 2D[d t + !(!y + V • b)]. 

Our new closed system of equations is badly nonlin- 
ear and coupled, but its linearisation can be immediately 
given in terms of an adjoint pair of Schrodinger equations 
with a potential ft, Indeed, id t tp = -DAip + ^tjj 

with a solution ip — p 1 / 2 exp(iS) and its complex adjoint 
makes the job, if we regard p together with S to remain 
in conformity with the previous notations. The choice 
of ^){x, 0) gives rise to a solvable Cauchy problem. This 
feature we shall exploit in below. Notice that, for time- 
indepedent il, the total energy J R s{\ — Q + typd 3 x of 
the diffusing ensemble is a conserved quantity. 

The general existence criterions for Markovian diffu- 
sion processes of that kind, were formulated in Ref. 
see also [plpl- 

Let us consider a simple one-dimensional example. 
In the absence of external forces, we solve the equa- 
tions (in space dimension one) dtp = — V(up) and 
dtv + (vV)v = +VQ, with an initial probability density 
Po(x) chosen in correspondence with the previous free 
Brownian motion example. We denote a 2 = ADto- Then, 
p(x,t) — [ 7r ( Q ,4 +4 £)2 t 2-)]i/2 exp[— Q 4 + £)2 t 2 ] and b(x , t) — 

v(x, t) + u(x, t) = 2D ai+4^if2 X are the pertinent solu- 
tions. Notice that u(x, 0) = —^§r- = b(x,0) amounts 
to v(x, 0) = 0, while in the previous free Brownian 
case the initial current velocity was equal to —D\7ln pq. 
This re-adjustment of the initial data can be interpreted 
in terms of the counterbalancing (recoil) phenomenon: 
the would-be initial Brownian ensemble current velocity 
vq = —uq is here completely saturated by the emerging 
forward drift bo — uq, see e.g. also ||. We deal also 
with a fictituous "repulsive" force, which corresponds to 
the compression (pressure upon) of the Brownian en- 
semble due to the counter-reaction of the surrounding 
medium. We can write things more explicitly. Namely, 
now: Q(x,t) = ^r^ifjif^ - 1) and the corre- 
sponding pressure term (VQ = ^VP) reads P(x, t) = 
— ^j^pgp p{x, t) giving a positive contribution +VQ to 



the local conservation law. The related Hamilton- Jacobi 
equation d t S + 5|V5| 2 = +Q is solved by S(x,t) — 

i^lwL-i — D arctani— 2 ^-). With the above form of 
Q(x,t) one can readily check that the Cameron-Martin- 
Girsanov constraint euqation for the forward drift of the 
Markovian diffusion process with a recoil is automatically 
valid for = \lnp+ S: 2Q = 2D[d t 4> + + V • 6)]. 

In anology with our free Brownian motion discussion, 

let us observe that presently < x 2 >= \ + 2D a '* ■ ^ 1S 
easy to demonstrate that the quadratic dependence on 
time persists for arbitrarily shaped initial choices of the 
probability distribution pq(x) > 0. That signalizes an 
anomalous behaviour (enhanced diffusion) of the perti- 
nent Markovian process when = i. e. il r — 2Q. 

We can evaluate the kinetic energy contribution 
SnP\dx — a 2( a 4 4^4 C 2 t 2) which in contrast to the Brow- 
nian case shows up a continual growth up to the terminal 
(asymptotic) value This value was in turn an initial 
kinetic contribution in the previous free Brownian expan- 
sion example. In contrast to that case, the total energy 
integral is now finite (finite energy diffusions of Ref. ) 
and reads J R {\v 2 — Q)pdx = ^ (it is a conservation 
law) . The asymptotic value of the current velocity v ~ f 
is twice larger than this appropriate for the Brownian 
motion, v ~ Jr. 

It is easy to produce an explicit solution to (7), (8) 
in case of f2(x) = \ r ) 2 x 2 — Dj, cf. pi, with exactly 
the same inital probability density po(x) as before. The 
forward drift of the corresponding diffusion-type process 
does not show up any obvious contribution from the har- 
monic Smoluchowski force. It is completely eliminated 
by the Brownian recoil scenario. One may check that 
b(x, 0) = -^r = "(a:, 0), while obviously b = ^ = - jx 
would hold true for all times, in case of the Smoluchowski 
diffusion process. Because of the harmonic attraction 
and suitable initial probability measure choice, we have 
here wiped out all previously discussed enhanced diffu- 
sion features. Now, the dispersion is attentuated and 
actually the non-dispersive diffusion-type process is re- 
alised: < x 2 > does not spread at all despite of the intrin- 
sically stochastic nature of the dynamics (finite-energy 
diffusions of Ref. || ) . 

It is clear that stationary processes are the same both 
in case of the standard Brownian motion and the Brow- 
nian motion with a recoil. The respective propagation 
scenarios substantially differ in the non-stationary case 
only. 
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